Self-energy corrections to anisotropic Fermi surfaces. 
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The electron-electron interactions affect the low-energy excitations of an electronic system and 
induce deformations of the Fermi surface. These effects are especially important in anisotropic 
materials with strong correlations, such as copper oxides superconductors or ruthenates. Here we 
analyze the deformations produced by electronic correlations in the Fermi surface of anisotropic 
two-dimensional systems, treating the regular and singular regions of the Fermi surface on the same 
footing. Simple analytical expressions are obtained for the corrections, based on local features of the 
Fermi surface. It is shown that, even for weak local interactions, the behavior of the self-energy is 
non trivial, showing a momentum dependence and a self-consistent interplay with the Fermi surface 
topology. Results are compared to experimental observations and to other theoretical results. 
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I. INTRODUCTION 

Anisotropic materials present different physics at dif- 
ferent energy scales, and their behavior or response to 
external probes is difficult to interpret. A big amount 
of experimental works have made possible to study the 
puzzling electronic properties of many anisotropic ma- 
terials which, in general, present potential technological 
applications. More theoretical effort is needed in order 
to understand the detailed experimental data which re- 
veal an unconventional behavior. In conventional met- 
als, the excitations that govern their low-temperature 
physics present well defined momenta lying at the three- 
dimensional Fermi surfaces. In the anisotropic materi- 
als, as layered transition metal oxides, unusual electronic 
properties appear and, under certain conditions, changes 
of the effective dimensionality occur. The electronic in- 
teraction effects are enhanced as the dimensionality de- 
creases and can change the fundamental properties of 
the material^. Therefore, due to both the anisotropy 
and the periodicity along the axis perpendicular to the 
planes, specific collective excitations appear absent in 
two-dimensional (2D) and three-dimensional (3D) elec- 
tron gases^. 

The high-temperature cuprate superconductors are 
among the most studied layered transition metals ox- 
ides, treated as 2D systems in many approaches, due 
to its strong anisotropy. In the hole doped cuprates 
the FS topology changes with doping from hole-like to 
electron-like^,. Recently a change in the sign of the 
Hall coefficient has been reported for heavily overdoped 
LaSrCuO^. The evolution of the FS in electron-doped 
copper oxide superconductors with doping has been re- 
ported by ARPES experiments to change from electron- 
pocket centered at the (tt, 0) point of the Brillouin zone at 
low doping to a hole- like FS centered at (tt, tt) at higher 
doping 7 . 

Other transition metal oxides as cobaltates or ruthen- 
ates are multiorbital systems and their FS present a com- 
plex topology with different sheets derived from the dif- 



ferent bands at the Fermi energy. A correlated 2D mate- 
rial particularly interesting is the S^RuCi, a ruthenate 
considered a model Fermi liquid (FL) system with impor- 
tant electronic correlations which have to be taken into 
account when interpreting photemission spectra^ to ob- 
tain a clear picture of the electronic properties, especially 
in the vicinity of the Fermi energy. In Sr2Ru04 the FS 
separates into three sheets a, (3, and 7, coming from the 
d X z, d yz and d xy orbitals. 

In the study of the electron-electron interactions in 
anisotropic metallic systems an open question is the de- 
formation of the Fermi surface induced by these inter- 
actions. The Fermi surface is one of the key features 
needed to understand the physical properties of a mate- 
rial and its shape provides important information. Re- 
cent improvements in experimental resolution have led to 
high precision measurements of the Fermi surface, and 
also to the determination of the many-body effects in the 
spectral function, as reported by angle resolved photoc- 
mission spectroscopy (ARPES) experiments 9 . However 
the interpretation of the data obtained by different ex- 
perimental techniques in anisotropic strongly correlated 
systems remains a complex taski^. 

The Fermi surface depends on the self-energy correc- 
tions to the quasiparticle energies, which, in turn, de- 
pend on the shape of the Fermi surface. Hence, there 
is an interplay between the self-energy corrections and 
the Fermi surface topology. For weak local interactions, 
the leading corrections to the Fermi surface (FS) arise 
from second order diagrams. The self energy, within this 
approximation, can show a significant momentum depen- 
dence when the initial FS is anisotropic and lies near hot 
spots, where the quasiparticles are strongly scattered^. 
This simultaneous calculation of the FS and the second 
order self-energy corrections is a formidable task. How- 
ever, the knowledge of the exact shape of the FS of a 
material is very important since it may affect the trans- 
port properties as well as the collective behavior, and 
have a valuable information from the point of view of 
theory in orther to find the appropriate model to study 
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the system. Many approaches have been used to study 
this problem like mean- fiel d 12 ' 13 , pertubation theory^, 
bosonization method o 15 ' 16 , or perturbative Renormaliza- 
tion Group calculation o 17 ' 18 ' 19 ' 20 ' 21 , and the cellular dy- 
namical mean-field theory (CDMFT), an extension of 
Dynamical Mean Field Theory^ 2 -, and many others. In 
spite of the great theoretical effort done in the last years, 
there is a need to develop alternative new methods in or- 
der to understand the origin of the electronic properties 
in materials with strong correlations. 

In this work, we calculate perturbative corrections and 
use Renormalization Group argument a 23 ' 24 in order to 
study analytically the qualitative corrections to the shape 
of the FS induced by the electron-electron interaction. 
This method allows us to classify the different features 
of the FS from the dependence of the self-energy correc- 
tions on the value of the high energy cutoff, A, defined at 
the beginning of the Renormalization process. As it will 
be shown later, one can also analyze the effects of varia- 
tions in the Fermi velocity and the curvature of the non 
interacting FS on the self-energy corrections. The calcu- 
lations do not depend on the microscopic model which 
gives rise to a particular Fermi surface, so that it can 
be useful in different situations. For concreteness we will 
consider the t — t' Hubbard model to study two dimen- 
sional Fermi surfaces of cuprates and an extension of it 
to study the case of Sr2Ru04. The paper is organized as 
follows. We define the model in Section II and describe 
the way the self-energy corrections are calculated. In Sec- 
tion III we present a detailed calculation of the changes 
expected for a regular FS, as well as for a FS showing 
singular points like Van Hove singularities, nesting or in- 
flexion points. We compare with results from ARPES 
experiments on anisotropic materials, mainly cuprate su- 
perconductors and Sr2Ru04. In the last section we high- 
light the most relevant aspects of our calculation, and 
some conclusions arc presented. 



II. THE METHOD 

The method of calculation of the self-energy correc- 
tions does not depend on the microscopic model used 
to obtain the electronic structure and the FS. In our 
scheme simple analytical expressions of the effects in- 
duced by the interactions are deduced from local fea- 
tures of the Fermi surface, and we are able to treat, on 
the same footing, the regular and singular regions of the 
FS. Therefore the method is particularly useful in corre- 
lated anisotropic materials which present exotic proper- 
ties and deviate from band structure calculations. The 
importance of considering correlation effects when inter- 
preting experimental data is already known and recently 
a great effort has been made in order to evaluate the self- 
energy from ARPES spectral. The evaluation of many 
body effects in these complex materials is far from triv- 
ial since the electron scattering presents a dependence on 
momentum and energy. We limit the study to the weak 



coupling regime, considering weak local interactions, con- 
sistent with the Hubbard model. 



A. The model. 

We consider the t — t' Hubbard model which is the 
simplest theoretical model which allows us to study dif- 
ferent correlated materials and describes the shape of the 
FS observed by ARPES in different materials as cuprates 
(see refj£ and references therein). Depending on the ratio 
t' jt and on the band-filling, different phases and stabil- 
ities appear, as found in early mean-field and Quantum 
Monte Carlo (QMC) studies of the model 2 ^. By changing 
the parameters a rich phase diagram, including antifer- 
romagnetic, ferromagnetic and superconducting phases, 
has been found for the 2D t — t' Hubbard model which 
describes many physical features of copper oxides and of 
Si^RuO^^l. For the cuprates, the most studied model is 
the Hubbard model on a square lattice and considering 
an effective single band. The Hamiltonian of the t — t' 
Hubbard model is: 



U J2 



(1) 



where c St i(c si ) are destruction (creation) operators for 

electrons of spin s on site i, n, lS = c' si c St i is the number 
operator, U is the on-site repulsion, and t,y = t are the 
nearest and tij = t' the next-nearest neighbors hopping 
amplitudes, respectively. The Fermi surfaces of the non 
interacting systems are defined by: 

£f = e(k) = 2t [cos^a) -I- cos(k y a)]+4t' cos(k x a) cos(k y a) 

(2) 

where a is the lattice constant. 
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FIG. 1: (Color online) Qualitative picture of the evolution 
of the FS with filling from almost isotropic to convex, going 
through a FS exhibiting inflexion points, and one with van 
Hove singularities (left panel, t' — — 0.3t). A region with 
almost perfect nesting is shown in the right panel {t! = 0.3t). 

Assuming that t < 0,t' > and \2t'\ < \t\, the Fermi 
surface is convex for -2t+W <e F <e = -8i'+16i' 3 /i 2 . 
For -8t' + 16t' 3 /< 2 < e F < -At' the Fermi surface shows 
eight inflexion points, which begin at k x = k y — kg = 
a -1 cos~ 1 (— 2t'/t) and move symmetrically around the 
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(±1,±1) directions, towards the center of the edges of 
the square Brillouin zone, (0,±7r), (±7r, 0). For e F = At' 
the Fermi surface passes through the saddle points (Van 
Hove singularities) located at these special points of the 
Brillouin zone. For At' < ep < — At, the Fermi surface 
is convex and hole like, centered at the corners of the 
Brillouin Zone, (±7r, ±7r). In Fig.pQ the variation of the 
FS shapes with doping is qualitatively shown. When only 
nearest-neighbor hopping is considered, t' = 0, the model 
has particle hole symmetry, and the Fermi surface shows 
perfect nesting for ep = 0. FS shapes similar to these 
shown in Fig.[T] have been experimentally observed by 
ARPES on different cuprate samples, at different doping 
levels. 



B. Self-energy analysis. 

We will analyze the interplay between the electron- 
electron interactions and the FS topology in the weak- 
coupling regime. The corrections to the non-interacting 
Fermi surface are given by the real part of the self- 
energy. For each filling n, a Fermi surface is defined. 
The electron-electron interaction leads to a self-energy 



S(k, to) = Re£(k, uj) + Tm£(k, uj) 



(3) 



which modifies the bare one-particle propagator 



Go(k,cj) 1 =u — e^+idsgne^ (where S — ► + ) to 



G(k,w) = 



1 



oj - [e(k) - e F ] - S(k,w) 



(4) 



and the FS of the interacting system is given by the u> = 
solution of the equation 



e F -e(k)-ReE(k,w) =0 



(5) 



As explained above there are many possible shapes of 
the FS which fit the experimental results from copper- 
oxide or ruthenate samples. The conventional perturba- 
tion theory fails on describing FS for which logarithmic 
divergences in the density of states (DOS) appear at cer- 
tain values of the parameters of Eq. {T]) . Then we proceed 
to calculate the self-energy by adopting a Renormaliza- 
tion Group strateg y 23 ! 24 . It is assumed that the effect of 
the high energy electron-hole pairs on the quasiparticles 
near the Fermi surface have been integrated out, lead- 
ing to a renormalization of the parameters t, t' and U of 
the hamiltonian. The possibility that other couplings are 
generated in the system is not allowed. Thus, the Hamil- 
tonian, Eq.JT]), describes low temperature processes be- 
low a high energy cutoff, A <C t,t'. For consistency, we 
consider the Hubbard interaction U < A, as well below 
the energy cutoff. Therefore the corrections to the quasi- 
particle energies are determined as a function of A, which 
defines an energy scale about the Fermi line which will 
contain the modes we are interested on (low energy ex- 
citations with | eg | < A), separated from the high energy 
excitations (with |eg| > A) which will be integrated out. 
We have to notice that we are restricting ourselves to a 
momentum-independent coupling U which corresponds 
to a local interaction in the real space. During the pro- 
cess we assume that the Fermi surface of the interacting 
system exists, and that this FS, dressed by the correc- 
tions due to the interactions, has the same topology as 
that of the non-interacting system. 

The two- loop self-energy shown in Fig.[2](b) can be 
computed from 



i£ 2 (k» = JdLo'Jd 2 qG a (k-q,LU-Lu')Il(q,Lo') 

(6) 

where the one-loop particle-hole polarizability, in 
terms of the one-particle propagator, reads 



where ReS(k,cj) is the real part of the self-energy. The 
diagrams that renormalize the one-particle Green func- 
tion up to second order in perturbation theory are de- 
picted in Figure [2]. The Hartree diagram, shown at the 
left of the figure, gives a contribution which is indepen- 
dent of momentum and energy, hence it cannot deform 
the FS. The two-loop diagram (right of Fig. [2]), modi- 
fies the FS through its k dependence and, in addition, 
it changes the quasiparticle-weight through its uj depen- 
dence. 




(a) 



(b) 



FIG. 2: Low order self-energy diagrams. Left: Hartree dia- 
gram. Right: two loop correction. 



ill(q,u) = -^ jduj'Jd 2 k@(A-\e^\)Q(A-\e^\)G (k,uj')G ^ + ^UJ + uj') (7) 

I 

The cutoff in energies A is used to implement the RG scheme^ 3 -: the virtual states in the loop of the diagram 
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shown in Fig.[2](b) have to be kept in the energy range 
determined by the cutoff. 



III. RESULTS 

After we compute the self-energy as explained above, 
we will analyze the deformation induced in the Fermi 
surface shape. We are interested in anisotropic two- 
dimensional FS similar to those measured for the 
cuprates which present regions with different scattering 
rates. For clarity, we consider first Fermi surfaces with 
hot spots, which are regions where the scattering of the 
quasiparticles is strongly enhanced. At these points, the 
scattering could be singular giving divergences of the sus- 
ceptibility. Separately we address the deformations of 
regular FS, curved surfaces which do not have singular- 
ities. These surfaces present a scattering rate relatively 
weak. We will show that once the Fermi velocity and 
curvature of the non-interacting FS are known, we can 
evaluate the corrections to the FS shape. Even consider- 
ing a weak local interaction, without momentum depen- 
dence, the effects are strongly dependent of the location 
at the FS. 



A. Self-energy corrections to the Fermi surface 
near hot spots. 

In anisotropic materials, the FS can present regions 
or special points which are called hot spots where the 
quasiparticles become strongly scattered and their behav- 
ior deviates from the conventional Landau Fermi liquid. 
The FS of layered transition oxides, as have been shown 
above, changes with doping adopting different shapes 
which lie close to Van Hove singularities, or present 
nested flat regions or inflexion points combined with reg- 
ular sectors. 

The effects of the Hubbard interaction have been stud- 
ied when the Fermi surface presents hot spots both, 
near the perfect nestin g 28 ' 29 i 30 or near a van Hove 
singularity^^ 3 ^ 3 ^! 3 !^. The curvature of the FS 
has important implications in the properties of the sys- 
tem and the inflexion points, which separate regions 
where the curvature has opposite signs, may induce 
anomalous effect o 39 ' 40 . The crucial role that the FS 
geometry plays on the unusual physics of 2D systems, 
makes it desirable a deeper insight in the interplay be- 
tween it and measurable parameters. The functional de- 
pendence of the self-energy on the cutoff is different in the 
vicinity of the hot spots than in regular zones of the FS. 
Near the hot spots to be considered here, the dispersion 
relation satisfies, near the Fermi level: 



±- 



n x 1 r. 
Vpk± 



Van Hove 
nesting 



(8) 



where eg = e(k) — e-p, k± is the momentum perpendic- 
ular to the FS relative to hp, k± = (k — k^)j_, Vp is the 
Fermi velocity at any particular point and m x ~ m y . 

Unlike the usual quadratic dependence expected in a 
Fermi liquid, the frequency dependence of the imaginary 
part of the self-energy in a nested region of the Fermi 
surface, or at Van Hove singularities is known to be lin- 



Im£ 2 (k, eg) 



(9) 



At the FS parts away from the hot spots, the leading 
contribution to the two loop self-energy, when the Fermi 
surface is near a Van Hove singularity, comes from dia- 
grams where the polarizability bubble, n(q, u>) expressed 
in Eq.((7j), involves transitions near the saddle point) 41 -. 

Only particles in the vicinity of hot spots on the FS are 
strongly scattered and present an anomalously large life- 
time, while away from the hot spots the single particle 
lifetime follows the Landau's Fermi liquid theory energy 
dependency. At some values of the band filling the FS 
is near a nesting situation, as shown in Fig.pQ, then the 
polarizability at low momenta is similar to that of a one 
dimensional Fermi liquid, due to the flat FS regions. The 
susceptibilities can be written as: 



n(q» 



m*|q|- 



Van Hove 
nesting 



(10) 



where m* is an average of the second derivative of the 
bands at the saddle point. Note that, in both cases, 
the density of states is proportional to the inverse bare 



bandwidth W 



t-\t'- 



The imaginary part of the second order self-energy near 
the regular regions of the Fermi surface can be written 



as 41 -: 



ImX^k, eg) ~ / did dqlmll(q 



(11) 



where <7 max ~ |A|/«f, and vf is the Fermi velocity in 
these regions. By combining Equations (flT)]) and (TlTj) . we 
find: 



Im£ 2 (k,eg) tx < £ 



e 3 J 2 Van Hove 



nesting 



(12) 



According to Eq. (|12)) the usual Fermi liquid result 
is recovered for the regular parts of the Fermi surface 
near almost nested regions. This result arises from the 
fact that the small momentum response of a quasi-one- 
dimensional metal does not differ qualitatively from that 
predicted by Landau's theory of a Fermi liquid while, 
close to the Van Hove singularities, the energy depen- 
dence of the ImX^k, eg) presents anomalous exponents. 
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The effects induced by inflexion points have been ad- 



dressed 



•39,40 



, where the instabilities of anisotropic 2D 



systems are analyzed. Near the inflexion points, the dis- 
persion relation can be expanded about the Fermi level 
and satisfies 



Up fell 
Vpku 



b 2 k\ 



inflexion point 



special inflexion point ^ ' 



where hi is the momentum parallel to the FS relative 
to kp, feii = (k — ki?)||, kj_ is the momentum perpen- 
dicular to the FS relative to fc^, fc_L = (k — k^ and 
b\ 2 are constants. The special inflexion points lie along 
a reflection symmetry axis of the BZ — (the k x = k y = ka 
points). We use the techniques previously developed in 
Ref. [ 39 ' 4 ° ] to obtain the second order self-energy near an 
inflexion point: 



Im£ 2 (k, eg) oc tJ 2 



(14) 



In the case of an special inflexion point, where the 
Fermi surface changes from convex to concave and a pair 
of inflexion points are generated for ep = eo and k = 
(ko,ko) defined earlier, the imaginary part of the self- 
energy behaves as Im£ 2 (k, eg) oc eV 4 . 

Once the imaginary part of the self-energy is known, 
we can obtain the real part of the self-energy from it by 
means of a Kramers Kronig transformation. Although 
Iml] 2 (k, eg) has been given for u) = eg, since we are 
in the weak coupling regime U ~ A <C cf, where cf is 
of the order of the non-interacting bandwidth, then the 
imaginary part of the self-energy associated to a state 
with energy eg is only significant in an energy range 
— A + eg < uo < A + eg. We assume that one can ap- 
proximate ImEg(w) in this range by an expansion on 
(uj — e^)/W, where W is an energy scale of the order 
of the bandwidth in the non interacting problem, and 
keep only the lowest order term. This approximation ne- 
glects contributions from a region of energies centered 
around \oj — eg| ~ A and of width SA, which is, at most, 
a fraction of A. The contribution of the Kramers Kronig 
transformation performed in this region of energies is, at 
most, of order ImSg(A) and does not modify the depen- 
dence of ReSg(eg) on the local properties of the Fermi 
surface. Therefore we can obtain the real part of the 
self-energy from the imaginary part by a Kramers-Kronig 
transformation, and restricting the frequency integral to 
the interval < to < A, we obtain: 



Re£ 2 (k, eg) cx - 5 2 |A| 




Van Hove 
nesting 



(15) 



where the negative sign is due to the fact that it is a 
second order contribution in perturbation theory, and g 
is a dimensionless coupling constant of order U /W. The 



sign is independent of the sign of U in Eq.([T]). In the 
regular parts of the Fermi surface, Eq. (fT2"|) leads to: 



Re£ 2 (k,eg) 



-fl 



2|A| 



3/2 



w 1 / 2 
y w 



Van Hove 
nesting 



(16) 



where the additional powers in W arise from the m* and 
vf factors in the susceptibility, expressed in Eq. (fl"0|) . 

In the limit A/W — > 0, the different dependence on A 
of the self-energy corrections at different regions of the 
Fermi surface is enough to give a qualitative description 
of the changes of the Fermi surface. For instance, when 
the non interacting Fermi surface is close to the saddle 
point, k = a _1 (±7r, 0), a _1 (0, ±7r), the self-energy cor- 
rection is negative and highest in this region. Note that 
the logarithmic divergences in Eq. (|15p are regularized by 
the temperature or elastic scattering. 

At band fillings where the FS lies close to the Van Hove 
singularities, most part of the low energy states close to 
the Fermi energy, are around the saddle points (0, ±7r) 
and (±7r, 0) (see Fig.([T])). Strong screening processes 
arise due to the big density of states at these points, 
and if the chemical potential of the system is kept fixed, 
i.e., the system is in contact with a charge reservoir, the 
number of particles varies and the Fermi energy tends to 
be pinned at the Van Hove singularities^. Then, in order 
to remove the Fermi surface from a Van Hove point or 
nesting situation, a large number of electrons must be 
added to the regular regions. When the points of the FS 
near these hot spots are at distance k from the hot spot, 
the change in the self-energy needed to shift the Fermi 
surface by an amount 8k is, using Eq. (|15[) . 



<5£ oc g 2 A^ 



(17) 



with additional logarithmic corrections near a Van Hove 
singularity. Near the regular regions of the Fermi surface, 
a shift in energy of order <5£ leads to a change in the 
momentum normal to the Fermi surface of magnitude 
(5fc r og ~ 5T,/vf- The area covered in this shift gives the 
number of electrons which are added to the system near 
the regular regions of the Fermi surface. We find 



Sn 



eg 



> femaxA 5k 



(18) 



where fc m ax ~ a 1 determines the size of the regular re- 
gions of the Fermi surface. The value of 8n diverges as the 
Fermi surface moves towards the hot spot, k — > 0. Hence, 
the number of electrons needed to shift the FS away from 
the hot spot also diverges. This result has been obtained 
from calculations at fixed chemical potential where 
the presence of a charge reservoir is considered, with reg- 
ular self-energy corrections. This situation has particular 
interest when studying the physics of high-X^ cuprates, 
where doping of the Cu0 2 layers and interactions with 
the rest of the perovskite structure are important. The 
pinning of the Fermi level to the Van Hove singularity has 
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been investigated in the 2D t — t' Hubbard model by RG 
techniques 43 , taking into account the formation of flat 
bands due to the renormalization of the electron spec- 
trum. The pinning of the Fermi level to the Van Hove 
singularities is found without making use of a reservoir, 
and the chemical potential of the system remains practi- 
cally constant in a range of dopings near the Van Hove 
filling. 



Self-energy corrections to regular Fermi 
surfaces. 



In this section we study a 2D system at a band filling 
which yields a curved FS, slightly anisotropic, in the ab- 
sence of singularities. Near the Fermi surface, by choos- 
ing an appropriate coordinate system, the electronic dis- 
persion can be approximated by: 



VFk± + (3k,, 



(19) 



where eg 



e(k) — ep, fc_L is the momentum per- 
pendicular to the FS relative to kp, k± = (k — kp)j_, 
fen is the momentum parallel to the FS relative to /cf, 

fe|| = (k — ki?)||, vf is the Fermi velocity at any particu- 
lar point i>f = nj^ • Ve(k), and (3 is related to the local 
curvature of the Fermi surface b = ny • ^V 2 e(k)^ ny, 

by (3 = bvp/2. This expansion implies, without assum- 
ing any rotational symmetry, a FS locally indistinguish- 
able from a circular one, where the energy Eq . (|19|) would 
correspond to a radius kp = rriFVF, where we have re- 
named (3 — \/2mp, being mp the effective mass. The 
Fermi velocity vp and the FS curvature 6, are functions 
of t, t 1 , ep and the position along the Fermi line. We cal- 
culate the second order diagram of Fig. ©, assuming that 
the main contribution to the self-energy arises from pro- 
cesses where the momentum transfer is small, forward 
scattering channel, or from processes which involve scat- 
tering from the region under consideration to the oppo- 
site part of the Fermi surface, i. e., backward scattering 
(in Appendix A we give the expressions of the polariz- 
ability for these two channels). This assumption can be 
justified by noting that we are considering a local Hub- 
bard interaction, which is momentum independent, so 
that the leading effects are associated to the structure of 
the density of states. The processes discussed here are 
those joining the regions which have the highest density 
of states. 

From Eq.© we can obtain the imaginary part of the 
self-energy, which describes the decay of quasiparticles in 
the region under consideration and that it is independent 
of the cutoff A. The contribution from forward scattering 
processes is: 

ImSi(k>) = (20) 

The quadratic dependence of energy is expected, and 
consistent with Landau's theory of a Fermi liquid. This 
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FIG. 3: (Color online) Real and imaginary parts of the self 
energy as a function of the energy uj: full red line — ReE, 
dashed blue line ImE 



contribution diverges as vf — > 0, that is, when the Fermi 
surface approaches a Van Hove singularity, or as |6| — > 
which signals the presence of an inflexion point or nest- 
ing. The contribution due to backward scattering is ex- 
actly the same as that from forward scattering, Eq. (|20|) . 
with the same numerical prefactors. 

Using a Kramers-Kronig transformation, and integrat- 
ing in the interval < uj < A, we obtain: 



ReS^k, uj) ■■ 



3 U 2 a 4 1 



64 y^n 3 v 2 



A 2 + 2Auj + 2uj 2 \og 



A 



(21) 

From the experimental point of view the determination 
of the scattering rate (ImS(k, uj)) presents particular in- 
terest and big effort has been devoted in order to obtain 
it: by ARPES because of the momentum and energy re- 
solved measurements ^ 25 ' 44 ^ 5 and recently by electrical 
transport experiments at microwave frequencies^. 

The extraction of the correlation functions from the 
experimental data is a complicated task and, although 
many theoretical approximations exist, the computation 
of correlation effects is also difficult. From ARPES 
results in underdoped and optimally doped cuprate 
samples 45 an anisotropic scattering rate around the 
Fermi surface have been found and the bare Fermi ve- 
locity has been directly obtained. By using a different 
methodology the real and imaginary parts of the self- 
energy has been obtained from photoemission data by a 
self-consistent procedure^ 5 -. 

In Fig. [3] we represent the self-energy as a function 
of the frequency according to results from Eqs. (|2QI2ip . 
where the linear (quadratic) behavior of the real (imag- 
inary) part of the self-energy at low frequencies, typi- 
cal of a Fermi liquid system, is recovered. We find a 
qualitative agreement with the low-energy part of the 
self-energy functions extracted self-consistently from the 
experiments in Ref. [ 25 ' 47 ] . It should be notice that we 
consider here the electron-electron scattering only. The 
impurity and electron-phonon scattering will, no doubt, 
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cause finite life-time and energy renormalization of the 
excitations but our main concern is the self-energy due 
to electro-electron correlation. The impurity scatter- 
ing term can be considered to be isotropic (from an 
isotropic distribution of static impurity scatterers) an it 
will give a constant term in Im£(k, ui). The electron- 
phonon self-energy can be assumed to be small at low 
temperature. Then, the assumption that the dominant 
scattering mechanism is the electron-electron interaction 
in the systems under study is not unreasonable 8 . The 
essential features of the electron-phonon coupling are ex- 
plained in Ref.[— ] both in the superconducting and nor- 
mal states. An effective ReE is extracted from ARPES 
measurements and by fitting and modelling of the data 
information about bosonic excitations. 

As a last result, from the real part of the self-energy 
Eq. ([2"Tj) . we can calculate the quasiparticle- weight 



(-it.it) 



(it, it) (-ji,h) 



(11,11) (-TT,7l) 



(It.rc.) 



z c 



<9Re£(k F ,w) 



dw 



(22) 



which for our case reads 



Zc = 

k F 



1 



64tt 3 



U 2 A 



(23) 



In the weak coupling regime, the quasiparticle-weight 
would be minimum either if the Fermi velocity becomes 
very small (Van Hove singularity, consistent with re- 
sults in Ref.f^]) or if the curvature of the Fermi surface 
changes sign (inflexion point). 

Finally, from Eq. (|2ip we obtain the expression which 
gives the zero frequency limit of the real part of the self- 
energy 



ReS 2 (k,w = 0) = - 



3 U 2 a 4 A : 
64 V2tt 3: 4[ 



(24) 



which renormalizes the FS according to Eq. ([5]) . 



C. Application to Fermi surfaces of copper oxides 
superconductors 

As stated in the Introduction, copper oxides super- 
conductors have a strongly anisotropic layered structure. 
Cuprates present anomalous properties in many physical 
aspects and are one of the main challenges to condensed 
matter physics. The peculiarities of the phase diagram 
have been addressed by different techniques and no con- 
sensus has been reached. Experimental results indicate 
a behavior far from the Fermi liquid, and the changes in- 
duced by doping on the ground state character of the nor- 
mal state add complexity to the problem. Both, the low 
dimensionality (CuOi planes) together with the strong 
electronic correlations, have to be taken into account to 
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FIG. 4: (Color online) Deformations induced by the interac- 
tions on the FS of the t — t' Hubbard model. The axes are 
labelled in units of a" 1 , where a is the lattice constant. Black 
line represents the unperturbed FS while the red line repre- 
sents the FS corrected by the interaction. For t' /t = —0.3 (a): 
high doping range and (b) close to half filling. For t'/t = +0.3 
(c) close to half filling. 



understand the low energy excitation spectrum of the 
cuprates. The effects of the strong correlations on the 
Fermi surface shapes of the cuprates is a hot issue. Re- 
cently Civelli et al2^ have addressed the problem by us- 
ing an extension of the dynamical mean-field theory, the 
cellular dynamical mean-field theory (CDMFT) which al- 
lows the study of A;— dependent properties. They study 
the 2D Hubbard model in the strongly correlated regime 
(U = 16t). A strong renormalization of the FS shape, 
due to interactions, is found together with a momentum 
space differentiation: appearance of hot and cold regions 
in the Brillouin zone. 



ReZ 




FIG. 5: (Color online) Real part of the self-energy for il — 
— 0.3t represented in the square Brilluoin zone. 



We study as well the two dimensional Hubbard model 
on a square lattice, and we will consider hopping am- 
plitudes tij to nearest-neighbors t and to next nearest 
neighbors t' . Since we are in the perturbative regime we 
consider local weak coupling instead of the strong cou- 
pling addressed in reference^. We adopt the hopping val- 
ues t = — 1 and t' — — 0.3t which mimic the hole-doped 
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cuprates (t/t' < 0), and we will consider two different velocity vp and the curvature b. For the dispersion re 



doping levels. In Eq. ([M|) we can see how the self-energy 
corrections due to electron-electron interactions depend 
on local features of the non- interacting FS, as the Fermi 



J 



lation given by Eq.([2), the expressions derived for the 
Fermi velocity and the curvature (taking for simplicity 
a = 1) at the Fermi level, are: 



t>i?(k) = \J (t + 2t' cos k y ) 2 sin 2 fc T + (t + 2t' cos k x ) 2 sin 2 fcj, 



(25) 



6(k): 



cos k y )[t cos k y sm 2 k x + 2t! cos 2 k y sin 2 k x + (t cos k x — t'(— 3 + cos 2k x )) sin 2 fcj,] . 
[(* + 2f cos k y ) 2 sia 2 k x + (t+ 2t' cos k x ) 2 sm 2 k v ] 3/2 



(t + 2t' cos k x )(t + 2t' cos k y 

[(t + 2f cosfc y ) 2 sin : 



r 



where k stands for the Fermi momentum k^?. These 
expressions illustrate the momentum dependence of the 
self-enegy corrections. At high doping, the FS presents 
an almost square shape with rounded corners. We find 
that the self-energy corrections are stronger for the re- 
gions with the smallest curvature at the diagonal parts 
of the BZ. This behavior, shown in Figure [1] (a) for 
e-p = —2.3 « 8t' , is similar to the result obtained by 
Freire et al in — for the renormalization of a flat FS by a 
two-loop field theory RG approach, where interactions in- 
duce a small curvature to the bare flat FS. They found as 
well that, the renormalized FS, becomes truncated due to 
the interactions, not found here. Next we will consider a 
lower doping level. By changing the filling, the FS shape 
varies, and close to half-filling, for ep = —0.9 = 3t', the 
FS has the form shown in Figure[?](b). The change in 
shape qualitatively agrees with the doping evolution of 
kp measured by ARPES on cuprates 4 ^, and the FS shape 
is similar to the FS reported in different experiments. At 
this doping, close to half-filling, the self-energy correc- 
tions enhance the hole-like curvature around (±7r, ±7r) 
and (±7r, =F7r), and flatten the FS close to the (±tt, 0) 
and (0, ±7r) points of the BZ as it is shown in Fig.[l](b). 
This result coincides with the renormalization found 
even though they are in the strongly correlated regime. 

Corrections found in Figure [3] can be understood by 
looking at Figure [5] , where the real part of the self-energy 
is depicted at u> = in the square BZ. At the central re- 
gion of the BZ near the T point, the main corrections 
occur around the diagonal of the BZ, as we have found 
in the high doping case, where the flat parts of the non- 
interacting FS become curved. At this point, the larger 
contribution to ReS(k, oj) will be due to the curvature of 
the Fermi surface, which is almost flat in the nodal region 
at this value of the band filling. In fact, it can be observed 
in Figure[5] , following the (— tt, — tt) — (tt, tt) diagonal line, 
that the minima of Re£(k,cj) coincide with the minima 
of the curvature, that corresponds to the maximum cor- 
rection to the non- interacting FS. 

Looking again at Figure [3] we see that, for higher fill- 
ings, the corrections are more pronounced where inflex- 



ion points start to appear in the Fermi surface (see Fig- 
ure Q]). The first of these inflexion points fco occurs in 
the diagonal of the BZ, in the nodal direction. Once 
this first inflexion point appears, if we increase the oc- 
cupation towards half-filling, new inflexion points merge 
in each FS and they distribute themselves symmetrically 
with respect to the diagonals of the BZ. This is why the 
divergence valley has this star-like structure. Finally, at 
lower doping levels, as that represented on the right panel 
of Figure [3], the Fermi line reaches the region closer to 
the border of the BZ, and the main corrections appear 
at the proximity of the antinodal points (0, ±7r), (±7r, 0). 
As can be observed in Figure[5], ReS(k,0) shows pro- 
nounced dips close to the saddle points, and therefore 
the FS is renormalized in this region. These minima are 
due to Van Hove singularities where the vp vanishes. 

This result agrees with one-loop functional RG calcu- 
lation of the self-energy in the weak coupling regime of 
the 2D t — t'- Hubbard model at Van Hove band fillings 49 
where vanishing of the quasiparticle- weight on approach- 
ing the antinodal points is found. Away from the Van 
Hove fillings a quasiparticle peak, with small spectral 
weight, emerges at (tt,0) and (0, n). 

The case of an electron doped system can be ana- 
lyzed by a particle-hole transformation of the Hamil- 
tonian which reverses the sign of t' , (t/t' > 0). For 
t! = +0.3t we find that, close to half-filling, the self- 
energy corrections are stronger in the proximity of the 
saddle points, where vp — > 0. As can be seen in Fig. [?] (c), 
the corrected FS is closer to a nesting situation than the 
bare FS. 

Our results, near half-filling, are in overall agreement 
with those of reference f 2 - 2 ], although we find that the self- 
energy corrections are stronger at the antinodal region in 
both, hole-like and electron-like, Fermi surfaces. 



D. Application to Fermi surface of Sr2Ru04 

Sr2Ru04 is a highly anisotropic layered compound, 
with an electrical anisotropy of about 400C&. It has 
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a strongly two-dimensional electronic structure and ex- 
hibits a good Fermi liquid behavior below 30K, as probed 
by bulk transport measurements^. The Fermi surface 
measured by ARPES matches the de Haas van Alphcn 
measurements and it can be well described by band struc- 
ture calculations. Therefore Sr2Ru04 is considered a 
correlated 2D material. As it occurs in the cuprates, 
the competition between superconducting and magnetic 
instabilities play an important role in the low energy 
physics of Sr2Ru04 which, with a critical temperature of 
about T c fts 1.5K, presents an unconventional supercon- 



ductivity with a p-wave and spin-triplet pairing^. This 
material has three relevant bands^, a, 0, and 7, of t2 g 
symmetry, formed from the 4<i-orbitals of the Ru 4+ ion. 
The {a, [3} bands are derived from {d xz ,d yz } orbitals 
and form two quasi-one-dimensional bands along the di- 
rections z and y respectively, that are weakly hybridized. 
The 7 band is derived from the d^-orbital and disperses 
into a real 2D band. In Figure[B] are depicted the elec- 
tronic structure, left panel and corresponding FS, right 
panel, calculated following the dispersion relation 



e 7 (k) = —2ti~ l {cos{k x a) + cos{k y a)) — At2- t cos{k x a) cos{k y a) — S{)~ f 
e a ,p(k) = — (*ic,/3 + ha,(})(cos(k x a) + cos(k y a)) 

±\J ((tia,p - t 2a ,p){cos(k x a) - cos(k y a))) 2 + 16t§^ sm(k x a) 2 sin(k y a) 2 - e 0a ,f3 (27) 



The tight-binding parameters are taken from 
references fl^]. This FS is in good agreement with that 
reported by ARPES experiments and obtained by band 
calculations. As before the self-energy corrections to 
the three sheets of the FS depend on the Fermi velocity 
and the curvature (not given here due to its size). 
In Figure [7] we can see the momentum dependence 
of the Fermi velocity for the three bands. It is easy 
to appreciate that the minima of the three plots are, 
besides in the T point, in the antinodal points, (±tt, 0) 
and (0, ±7r), and in the diagonal of the BZ, (±tt, ±tt) 
and (±7r, =F7r). 



eo —n 7i 




FIG. 6: (Color online) Left: Band structure of S^RuCU 
obtained from Eq. (|27|l using the parameter values given in 
Ref.—. The blue sheet corresponds to the Q-band, the red 
sheet is the /3-band and the green sheet correspond to the 
7-band. Right: Corresponding Fermi surface. 

The momentum dependence of the real part of the self- 
energy for the three bands is shown in Figure [5]. In the 
left hand side we can see the correction to the a-sheet, 
that is significant only either when the Fermi line lies 
near the T-point or in the proximity of the corners of the 
BZ, (±7r, ±7r) points. In a recent experimental work 8 -, the 
form of IrnS(uj) extracted for the bulk a-band from the 
ARPES spectra, is found to be consistent with a Fermi 



liquid, and that the quasiparticles residing in the surface 
layer a-band, show similar many-body interactions. (No- 
tice that the analysis and fitting procedure over raw data 
of this band made in Ref.f 8 -] leads to self-energy curves 
of the form of Fig.[3]). In the central figure, the contri- 
bution corresponding to the /3-band is shown. Here we 
can see that the correction is maximum (most negative) 
in the zones between the diagonals, due to the nearly flat 
regions of the Fermi line corresponding to the /3-band. 

Finally, the corrections to the 7-band are shown in the 
right hand side of Figure [5] (notice the different scales in 
the three graphs). In this case the corrections are max- 
ima near the antinodal points (±7r,0) and (0,±7r), due 
to the proximity of this band to a Van Hove point. In 
reference^ it is pointed out that calculated 7-band prop- 
erties depend very sensitively on how close it approaches 
the Van Hove points (7r,0), (0,7r). 

The main corrections occur at the 7-band as it is 
shown in Figure [9] where the bare and renormalized FS 
are depicted. Here non-interacting FS (full blue line) is 
changed, due to electron-electron correlations, to the in- 
teracting FS (dashed red line). This result agrees with 
photoemission measurements which indicate that the 7- 
band has much stronger interactions and plays a domi- 
nant role at low temperature^. 

The corrections in the a and (3 bands are due to cur- 
vature effects, while the corrections to the 7 band are 
due to Fermi velocity effects, because of the proximity of 
this band to the saddle points, as it has been pointed out 
above. Then, in agreement with ARPES measurements 8 
and other theoretical results^, the main corrections in- 
duced in the FS of Sr 2 RuC>4 occur at the 7-band cor- 
responding sheet. The importance of the FS geometry 
has been analyzed in Ref. f^] where the non-analytic cor- 
rections to the specific heat and susceptibility of a 2D 
Fermi liquid have been considered and the results ap- 
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FIG. 7: (Color online) Fermi velocity of the non-interacting bands: Left graph corresponds to the a-band, center graph 
corresponds to /9-band and right graph corresponds to 7-band of S^RuCU. 
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FIG. 8: (Color online) Real part of the self-energy corrections for the three bands of Sr2Ru04 in the first square BZ: Left graph 
corresponds to the a-band, middle to the /3-band, and right graph corresponds to 7-band. We have used the parameter values 
U — 0.01 and A = 1. Notice the different scales in the three graphs. 



plied to Sr 2 Ru04. Both, the dependence of the 7 band 
properties on how close the band approach the Van Hove 
points and the dominant interaction in the 7 band are 
found in Ref.f^], because it has the highest density of 
states at the Fermi level, as well as the largest mass and 
susceptibility enhancements^. 

Similarly, the importance of the band structure prop- 
erties of these materials can be seen in the context of 
multi-layer ruthenates, for which the proximity of their 
Fermi surface to a Van Hove singularity can give rise to 
a quantum critical end-point in the magnetic phase di- 
agram, as it is found, within a mean-field analysis, in 
Ref.[S£]. 

IV. CONCLUSIONS. 

The unconventional physics shown by anisotropic ma- 
terials, as high-T c superconductors and ruthenates, poses 
big difficulties when a theoretical model has to be chosen. 
On the other hand, although it is generally accepted that 
the many-body interactions of electrons play a key role 
in the underlying physics of these compounds and may 
be related to the occurrence of superconductivity in the 
cuprates, a consensus has not been reached about the 
origin of important features widely observed by differ- 
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FIG. 9: (Color online) Bare (continue blue) and interacting 
(dashed red) Fermi surface of Sr2Ru04. 

ent experimental techniques. The anisotropy momentum 
space shown by many electronic properties of the planes 
adds complexity to the possible theories. 

In the cuprates, the pseudogap phase, metallic but 



11 



with a broken Fermi surface (segments known as Fermi 
arcs^i) is an example of the remarkable momentum de- 
pendence of the interactions. Furthermore the band 
renormalization observed by ARPES in different families 
of high T c superconductors, known as kink in the dis- 
persion, which indicates a strong coupling to a collective 
mode (see Ref.^] and references therein), shows different 
energy scales and temperature dependence at the nodal 
and antinodal regions of the BZ, suggesting two different 
kinks of different origin, which nature is under debate: 
both phonon and magnetic mode have been suggested as 
possible causes. 

More recently a high energy anomaly in the spectral 
function has been observed in three different families 
of high-T c superconductors, and apparently in several 
ruthenate compounds^. This anomaly indicates that the 
quasi-particles at cf are dressed not only by the interac- 
tions with bosons at low energy, but also by interactions 
at higher energies. This peculiar high energy behavior 
together with the unconventional low energy properties, 
poses a challenge on theoretical models. 

The knowledge of the dressed FS is crucial to un- 
derstand the behavior of correlated materials, especially 
when an effective model is needed to explain the uncon- 
ventional physics. We have presented here a simplified 
way of taking into account the self-energy corrections to 
the Fermi surface. We have made use of the different 
dependence of the self-energy on the high energy cutoff 
in order to analyze the main features of the changes of 
the FS. The analysis presented here is valid only at weak 
coupling, and we do not consider corrections to the inter- 
actions or to the wave-function renormalization. On the 
other hand, the expressions obtained are analytical and 
related to the local features of the non-interacting FS in 
a simple way, so that they can be readily used to get an 
estimate of the corrections expected. 

The results suggest that the main self-energy correc- 
tions, which are always negative in our scheme, peak 
when the FS is close to the (±7r, 0), (0, ±7r) points in the 
Brillouin zone. If these contributions are cast as correc- 
tions to the hopping elements of the initial hamiltonian, 
we find that the nearest neighbor hopping, t, is weakly 
changed (as it does not contribute to the band dispersion 
in these regions). The next nearest neighbor hopping, t' 
which shifts the bands by —At' in this region, acquires a 
negative correction. This implies that the absolute value 
of if grows when if > 0, or decreases, when if < 0, in rea- 



sonable agreement with the results reported inS 2 -. Note 
that the tendency observed in our calculation towards 
the formation of flat regions near these points, when an- 
alyzed in higher order perturbation theory, will lead to 
stronger corrections. Our results also confirm the pin- 
ning of the FS near saddle points, due to the interactions. 
The analysis presented here is consistent with the mea- 
sured Fermi surfaces of the cuprates*^ and qualitatively 
agree with the doping evolution reported by ARPES^^. 
The self-energy corrections found for the FS of Sr 2 Ru04, 
which mainly renormalize the 7 sheet, are as well in quali- 
tative agreement with ARPES measurements^ and previ- 
ous calculations^. The broad spectrum of experimental 
data available at this moment makes comparison between 
results from different techniques one of the most efficient 
methods to obtain information about response and cor- 
relation functions of unconventional materials. To get 
an estimation of the self-energy corrections to the Fermi 
surface in a simple way, independent of the model, as the 
one here proposed, is helpful in order to gain insight into 
many low-energy physics aspects. 
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APPENDIX A: POLARIZ ABILITY 

In this Appendix the particle-hole polarizability for 
the forward and backward scattering channels are given. 
These polarizabilities are calculated in oreder to ob- 
tain the self-energy Eq. ([2H)l . The imaginary part of 
II(q, u>) for the forward channel is, using Eq.([7]) and the 
parametrization from Eq. (ll9|) : 



ImIT 



u 2 



2ft 



16-7T ftVF 



M| 



(Al) 



where q is a small vector that connects two pieces close 
together in the FS and MS = v F q± + §/3qjj. The argu- 
ment of the square root has to be positive, what gives an 
extra condition, usgn(ft) > sign(/3). 

Similarly, for the backward scattering we obtain: 



Imn B (Q + q, lu) = - 



U 2 1 

167T ftv-p 



^2/?(M+ Mf)-^2/3(Mg-H ) if M < M§sgn(/3) 
sgn(ft) ^2/3(Msgn(/?)+Mf) 



M > |Mf I 



(A2) 



where q is the deviation of the wave-vector from the vector Q that connects the region studied with the op- 
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posite part of the FS, and M? = vpq_L — |/3qii- 

For q = we have that ImII B (Q + q, lu) ~ -y/M, which 
agrees with the results obtained in Ref . , where a spin- 
fluctuation model for the Q = 2/cf instability is studied. 



For small w and fixed q we have, expanding Eq. (|A2|) up to 
first order in M/M B (for the case \iu\ < M B sgn(/3)), that 

ImII B (Q + q, u) ~ |cj|, as expected for a Fermi-liquid^. 
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